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THE DIAGONAL TWO-POINT CORRELATIONS OE THE ISING MODEL ON THE 
ANISOTROPIC TRIANGULAR LATTICE AND GARNIER SYSTEMS 

N.S. WITTE 


Abstract. The diagonal spin-spin correlations of the Ising model on a triangular lattice with 

general couplings in the three directions are evaluated in terms of a solution to a three-variable extension 
of the sixth Painleve system, namely a Gamier system. This identification, which is accomplished using the 
theory of bi-orthogonal polynomials on the unit circle with regular semi-classical weights, has an additional 
consequence whereby the correlations are characterised by a simple system of coupled, nonlinear recurrence 
relations in the spin separation N 6 Z>o. These later recurrence relations are an example of the discrete 
Gamier equations which, in turn, are extensions to the "discrete Painleve V" system. 


In this study we investigate the spin-half Ising model on the anisotropic, homogeneous triangular 
lattice and in particular the diagonal spin-spin correlations of the model. We identify a particular 
three-variable Garnier system in Corollary [i] as the integrable system lying behind this model and 
thus find an extension of the evaluation by Jimbo and Miwa in 1980 m of the same quantity for the 
square lattice case in terms of a t function of the sixth Painleve system. This extension is quite literal 
in the sense that the anisotropic triangular model includes the diagonal and row or column correlation 
functions of the rectangular lattice as special cases. To our knowledge this is the first appearance 
of a Garnier system in a statistical mechanical model even though such systems have arisen in other 
mathematical physics contexts: see fl 22 l for the ADM approach to 2-|-l-dimensional gravity; m for the 
symmetry reductions of the self-dual Yang-Mills equations in 2-1-2 dimensions; j 3 ^ for a generalised 
Henon-Heiles system; and iizi for finite gap solutions to the Garnier system and the g-dimensional 
anisotropic harmonic oscillator in a radial quartic potential. 

The first identification of a Painleve system for the Ising model correlations was for the scaled 
and critical correlation^ for the square lattice and in this case was made to a special type of the 
Painleve III system by Wu, McCoy, Tracy and Barouch in The extension to the lattice model for 
all T and finite N was found by Jimbo and Miwa m- This latter result for the spin-spin correlation 
{<^ 0 , 0 <^n,n) can be summarised as 
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Here t = k~^ and k = smh2i<Ci srnh 2 i<C 2 where Ki,K 2 are the couplings between sites along the hor¬ 
izontal and vertical edges respectively. This equation is solved subject to the boundary condition as 
t —)• 0, i.e. in the T <Tq regime 
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Where the limits N —>■ co and T —>■ Tc are carried out together in a certain combination. 













Alternative derivations of this result can be found in 13 ) and l 7 ). 

Equation Il2)l is a specialisation of the general n-form of Painleve VI which is constructed from the 
Hamiltonian system for Painleve VI llldl . Ili^ 

dq _dH dp _ dH 

dt ~ It ~ 

where the non-autonomous Hamiltonian is given by 

( 5 ) - 1)^^ = '?('? - 1)(^? - - K(‘? - 1)(^? - f) + - 0 + - 1)]P + a 2 (ai + ‘^ 2 ){q - 0- 

Here the D4 root system constraint ocq + cci + 2a.2 + ^3 + ^4 = 1 applies. The general n-form of the 
sixth Painleve equation is then 

(6) h'{t) [f(l - t)h"{t)f + [h'{t){2hit) - {2t - l)h\t)) + hb2hbif = Ylih'it) + bl), 

k=l 

and h{t) is related to the Hamiltonian by 

(7) h{t) = t{t - + e2{bi,b5,b4)t - \e2{bi,b2,b2„bi), 

where e2(') is the second elementary symmetric pol5momial of its arguments. 

In the present work we require the theory for the multi-variate extension of the Painleve VI system. 
Let us place M -|-1 > 3 independent variables in canonical position (see i-e. taking them at 

distinct points 

(8) to = 0, tl/ • • • / tM-2/ Im-I = 1/ fM = 00, 


each of which is associated with the parameters po/Pi/• • ■/i0M-2/|0M-i = P>Pm = Poo respectively. 
The Gamier system Qm-2 of typo L( 1 ^;M — 2 ) is labelled by an abbreviated Riemarm-Papperitz-like 
symbol (see §4 of IH) for the definition) in the following manner 

I to = 0 tl • • • tM-2 tM-1 =1 00 I 

(9) 


eo = n-po 01 =-pi ••• 0m-2 = -Pm- 2 0=-p 0OO = 2n-Fl - E^oV; 


The dynamics of the Gamier system is governed by the Hamiltonian system {qj,Pj;Kj,tj}^^^ with 
co-ordinate qr and momenta pr and with the Hamiltonian ISllHISD 
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The monodromy exponents are given by 9 j = —pj for ; = 1 ,..., M — 1 , 0 o = w — po/ arid 600 = 
2 n -|- 1 — Ej^o^P; with the constant k = —n(l -|- mo). The polynomial W, the denominator spectral 
polynomial, is given by 


(11) 


M -2 M -1 

W(z) = z(z - 1 ) n (Z - tj) = z ^ (-) 
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where the elementary symmetric functions of the singularity positions are denoted £;, Z = 0 ,..., M — 1 
and in particular eo = 1 / = 0 and Cm-i = Hj^i^ tj. The polynomial 2 V, the numerator spectral 

polynomial, is 
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(12) 2V(z) =z(z-l) (z-ty) 

7=1 

where the last relation defines the coefficients m;, I = 0 ,..., M — 1 and we observe that mq = po + 
P + Ej^T^P; and ?km-i = Po^M-i- The remaining polynomial, termed the spectral coefficient of 
degree M — 2 , is 
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This Hamiltonian system llTot is not polynomial in the canonical variables qr and the dynamical 
equations for qr do not possess the Painleve property in tj, however using the well-known canonical 
transformation to the new Hamiltonian system T-Lm-i = ESI/ pql 


(14) 

(15) 

(16) 


^ ^ r=l — 1) ((Jr — f;)©n(<?r) 


both these deficiencies can be removed. The sixth Painleve system is the M = 3 case of the above. 

The n-form differential equation ||2)l is quite useful in certain types of analysis of the correlations 
however additional insight can be gained from a recurrence sysfem in N which generafes the correla¬ 
tion at N -|-1 from ones at N or earlier, see for example the use made of these in |j^. In Propositions [1] 
and [2] a system of fwo coupled, first order non-linear difference equations for the diagonal correlations 
of the square lattice Ising model is given as (a = k~^) 
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although different but entirely equivalent systems have been given in fS) and The above system 
is a specialisation of the general "discrete Painleve V" system lIToB , [27! or the system in Sakai's 
classification m 


(19) 


r ^ [(v„ + n-t-po{t + l) - {pt + pi)t][wn + n - t - po{t + 1) - pt - pit] 
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(20) cvn + (^n-1 + (2w — l)f — 2 — 2po(f +1) — 2pt — pi{t + l) = {n — po)-z -- -|- (n -|-1 -h "^0)77-7- 
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All of the features that apply in the diagonal correlations of the square lattice case as a consequence 
of the identification with an integrable system, also apply in the case of the triangular lattice but only 
in a slightly more complicated way. Thus the main results of our study are the systems of non-linear 
recurrence relations for the diagonal correlations of the triangular Ising model in Corollary and 
and by implication for the column/row correlations of the anisotropic square lattice Ising model in 
Corollary and Corollary [7] We will not undertake the task of writing down nor analysing here the 
coupled, partial differential equations that follow from i fTot and which directly generalise ([Hi as this 
would divert us from our primary goal. 

The Ising model on the triangular lattice has attracted some considerable interest, partly because 
it exhibits the phenomena of frustration in the anti-ferromagnetic phase and partly because it led to 
the simplest form for the diagonal correlations on the square lattice. The partition function has been 
evaluated through all the methods known to work in treating the Ising model: Husimi and Syozi lITSi . 
[361 used a simplified version of Onsager's algebraic method on the isotropic lattice; Warmier 
used the method of Kaufman and Onsager for the isotropic lattice whereas Temperley [37! employed 
this approach on the anisotropic lattice; Newell first employed Kaufman's method in [23] and then 
Kaufman and Onsager's method on the anisotropic lattice in I26I : Potts employed the combinatorial 
method of Kac and Ward on the anisotropic model in |3o| . In all of these works the isotropic lattice 
was observed to possess a Curie point singularity in the ferromagnetic case, just like that for the 
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square lattice, but no singularity was found in the antiferromagnetic case. Subsequently Stephenson 
produced a series of works I-IV, [321 [^[3^[33] , sfarfing with the Pfaffian representation of fhe partifion 
function introduced by Kastel5m and utilised by Montroll, Potts and Ward to compute the square 
lattice row/column correlations [2^ . It is Stephenson's first and fourth works that we will primarily 
draw upon here. In his first work the two-spin correlations are evaluated as a Toeplitz determinant; 
in the second some particular four-spin correlafions are evaluated; in the third the isotropic anti¬ 
ferromagnetic lattice is studied using the two and four-spin correlafions and the asymptotics foimd for 
fhese at large site separations; and in the fourth we find the most detailed treatment of the anisotropic 
lattices in both the ferromagnetic and anti-ferromagnetic phases. Subsequent to these earlier studies a 
number of modern re-derivafions of fhe parfition function (and extensions thereof) [2^ , |j^ have been 
made by exploding fhe free-fermionic characfer of the model using Grassmann variable techniques. 

The triangular Ising model is known to be in the Villain-Stephenson imiversality class with critical 
exponents a = 0 , (2 —a)/v = 1 ,v = 1,/3/v = 1 / 8 , S = 7,7/v = 7/Al,S./v = 7 / 4 , see for example 
M . As menfioned above because the triangular Ising model in the anti-ferromagnetic regime is a 
prototype for geometrical frustration an extensive modern literature has emerged which we will only 
briefly touch upon. For example there are the studies M' liSl,l 4 D of fhe Ising model on a triangular 
Kagome lattice which explore this issue. Another feature of this effect is that the ground state of 
the model is highly degenerate, ^configurations = and has residual entropy at T = 0 . A 

consequence of this is that the boundary conditions do affect this entropy as #sites —?■ 00 and some 
computations 1)23), El have clearly illustrated this. Perhaps more relevant to this work is the study | [^ 5 ] 
of the pair (both diagonal and off-diagonal) correlations of fhe isofropic triangular Ising model in the 
anti-ferromagnetic phase at T = 0 and T > 0 through a number of approaches: fhe exacf resulfs for 
N < 20 numerically evaluated, asymptotic approximations as N —)■ 00 and Monte-Carlo simulations. 
Here it should be noted that in the isotropic anti-ferromagnetic case the model has no long-range order 
for T > 0 , i.e. the critical point is at T = 0 . What is observed here is the lack of exact results even in 
this highly specialised case and for T > 0 the existing approximations are of limited accuracy. 


1. Correlations along the diagonal of the anisotropic triangular Ising model. The Ising model is 
a system of spins cy G {—1,1} locafed at site r = {i,j) on a triangular lattice of dimension (2L -|- 
1 ) X ( 2 L -l- 1 ), or equivalenfly with the sites on a rectangular lattice having in additional to nearest- 
neighbour couplings along the x and y axes a third next-nearest-neighbour coupling along the up-right 
diagonals. Our conventions for fhe labelling of sites and of fhe coupling constanfs i = 1 , 2,3 are 
displayed in Fig. [ij wifh fhe origin af the centre of the lattice. 
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Figure i. Conventions for the co-ordinate system of the spin sites {i,j) and couplings 
Ki,K2,Kj for the homogeneous, anisotropic triangular lattice Ising model. 


The probability density function for a configuration is the product of Boltzmarm weights 

along the three axes 


(21) P[{t^y}J;y=_L] = 


Z 2 L +1 


exp 


L L-1 L L-1 

E E 

j——Li——L i——Lj——L 


L-1 L-1 

‘^3 E E 

i=-Lj=-L 


where Z2L+1 is the partition function. Averages such as the order parameter or average magnetisation 
are defined by 

( 22 ) {(Tofl) = LL 

however our interest will lie in the spin-spin correlation function 

( 23 ) {cro,OCrN,N) = E ‘^0,0>^N,N'P[{cri,j}lj^_L]- 

In all the averages we will take the thermodynamic limit, so for example limL_s.c»(Go,o) = (c^)- Our 
independent variables will either be the ordered triple (zi,Z2,Z3) or (ui,U2, U3) which are related to 
the coupling constants X„ i = 1,2 ,3 by 

1 — u I 0f 2 j 

(24) Zi = tanhfC, = --, u, = 

where the physical variables Z; G [— 1 , 1 ], u, 6 [ 0 ,00) for i = 1 , 2 ,^ In our characterisation of these 
correlations as a classical solution to the Gamier equations we will see that they are meaningful for 
complex values of the variables Ki, K2, X3 G C U {00} modulo certain restrictions and in fact this is the 
natural setting for these systems. 

From the symmetries of the correlations, see Eq. (1.6) of [35!, we have the Toeplitz matrix element 


(25) Wn{K^,K2,K^) = Wn{-K^,-K2,K^), 

(26) ={-lY+^Wn{-Kr,K2,-K^), 

(27) ={-\Y+^Wn{Kr,-K2,-K2,), 


^This notation is the reverse of that adopted by Stephenson whereas our conventions conform to common usage, see [23. 
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from which Stephenson inferred a classification with two classes: Class A which can be transformed to 
a completely ferromagnetic lattice Ki,K2, K3 > 0 and class B which can be transformed to a completely 
antiferromagnetic lattice K^, K2, < 0 . In the class A system or ferromagnetic system there is a critical 

point, the Curie point Tq, given by 

(28) 1 +Z1Z2Z3 = Zi+Z2+Z3+Z1Z2+Z1Z3+ Z2Z3, ViV2 + ViV3 + V2V3 = l, 

and exhibits a phase transition with a low temperature, 0 < T < Tq , ordered phase where the spins 
are aligned (a) > 0 and a high temperature, Tq < T < 00, disordered phase with no alignment 
(cr) = 0 . By contrast, in the antiferromagnetic class or class B systems there are two critical points, the 
Neel point Tjv 

(29) 1 + Z1Z2Z3 = -Zi - Z2 + 23 + Z1Z2 - Z1Z3 - Z2Z3, ViV2-ViV3-V2V3 = l, 


and the disorder point To > Tjv 

(30) Z3 + ZiZ2=0, -1^11^2 + 5^1^3 + 5^21^3 = 1- 


In the low-temperature regime 0 < T < the system exhibits antiferromagnetic long-range-order 
along the two lattice axes with the strongest K values and ferromagnetic order along the third. There is 
now an intermediate regime Tjv < T < Tp in which there is antiferromagnetic short-range-order along 
the two axes with the largest absolute values of K and ferromagnetic short-range-order along the third. 
In the high temperature regime Td < T < 00 there is exponential decay of the pair correlations along 
all the axes. There is a clear exposition of the 3-D phase diagram in Igj with respect to the co-ordinates 
{K-[, K2, K3) G IR^, where it is clear there exist rays emanating from the origin with no critical point for 
T > 0 , i.e. the system is disordered down to T = 0 such as in the isotropic anti-ferromagnetic case 
K^=K2=K3< 0. 

The spin-spin correlations along the diagonal have the Toeplitz determinant form given by Eqs. 
(6.10) and (6.12) derived by Stephenson [ j^ , or Eq. (1.4) of m using the methods of Montroll, Potts 
and Ward |g] 

(31) {<^o,oO'n,n)'^ = det[z(;y_jt]j^j;=o,...,N-i/ 


where the Toeplitz matrix element is the trigonometric integral 

^ de acos{ne)-bcos{{n-l)e)-ccos{{n + l)e) 

J-n 27 t _|_ ^2 _|_ (,2 _ 2 a{h + c) cos{ 6 ) + 2 bccos{ 29 )' 

with parameters a, b, c depending in the three couplings zi, Z2, Z3 

(33) fl = 2z 3(1-|-Zj)(l-I-Z2)+4 ziZ2(1+Z3), b = zlc = zl{l - zl){l - zl). 


Eor our purposes we note the essential feature that these matrix elements are the Eourier coefficients 
of the weight 


(34) 


00 

= E 

n——cio 


ya-b^-'^-cU 


The above general result subsumes a number of special cases which are of interest in their own 
right: 

(1) Diagonal correlations on the square lattice K3 = 0 , 

(35) (c^o,oc^n,n)'^ = 


(2) Row correlations on the square lattice Ki = 0 , 

(c^O,0<2‘N,o)'^ = {‘^0,0O'N,n)^ 


(36) 








(3) Column correlations on the square lattice K2 = 0 , 

(37) = {o-0,0I^N,n)^ 


K 3 h^K 2 


(4) Curie Point T = Tq, 

(5) Neel Point T = Tjq, 

(6) Disorder Point T = Tp. 

These correlations are governed by a Gamier system because of the simple observation: the weight 
function is a regular semi-classical weight, which is made precise below, and are thus these Toeplitz 
determinants are known to be T-functions for such equations IZiiii- In fact we are going to heavily 
employ results from the theory of Gamier system specialised to the classical cas^ in the particular 
setting of bi-orthogonal polynomials and their associated fimctions defined on the unit circle with 
regular semi-classical weights. The general theory of such systems has already been formulated for 
systems with an arbitrary number of singularities, following earlier work on orthogonal polynomials 
on the real line with semi-classical weights in 


Corollary 1 (|[^,[|4T|). The Toeplitz determinants {(cro,o^N,N)^}N=0' |j^, are classical 

T-functions of a Gamier system of type L(l^; 3 ) in three variables formed from the ratios of any three of 
the fourth (assumed non-zero) and ivhose abbreviated Papperitz-like symbol is 

(38) 


I 0 ^2 ^3 U 

\ eo = N 01 = -1/2 02 = -1/2 03 = 1/2 04 = 1/2 

where the formal monodromy exponent 6 j corresponds to the singularity 


00 

0OO = N 


Proof The logarithmic derivative of the weight fimction 

1 


has the form 


" »(C) ^ 


* 5 ’’ a>(C) dC'’” W(f)' 

where W(^), 2 y(^) are irreducible pol5momials in It is clear that 

(40) W(0 = (^ - ^l)(^ - ^2)(^ - ^3)(^ - ^4) = + £2^' - £3^ + £4, 


(41) 


2V{0 = WJ2 


Pi 






= -h m2C - m3. 


where the weight data is represented in the following way when the four singularities j = 1,2, 3,4 
are in generic positions and one really has the case of five finite singularities M = 5 as the one at 
the origin is always present (as is also the one at infinity under generic circumstances). Here V{z), 
W{z) are irreducible polynomials that must satisfy the following generic conditions of the regular 
semi-classical class: 


(i) deg (W) = M > 2 , 

(ii) deg (P) < deg (W), 

(hi) the M zeros of W(z), {zq, .. .,Za4-i} are pair-wise distinct, and 
(iv) the residues pj = 2 V{zj)/W'{zj) i Z>o. 

The above generic case defines a system of bi-orthogonal polynomials and their associated functions on 
the unit circle, and these satisfy the same second-order linear differential equation which is represented 
by the following abbreviated Riemann-Papperitz-like symbol 


0 Cl C2 C3 C4 00 I 

00 = n 01 = -pi 02 = -pi P3 = -P3 P4 = -P4 0OO = 2 n -h 1 - Poo J 


^The term classical used in Painleve theory is distinct from that employed in orthogonal polynomial theory and refers to 
the fact that the parameters of the Gamier equation are located on a Weyl chamber wall, i.e. a particular integrality condition 
applies to the parameters and the monodromy matrices are either lower/upper triangular or trivial. 
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Here £4 7^ 0 and ( [40) and ( [41) still apply. In our example M = 5 , with four finite, non-zero regular 
singularities 


0 Cl C3 C4 00 1 

n Pi = 1/2 p2 = 1/2 p3 = -1/2 p4 = -1/2 n j ' 


and thus one of the conditions on the weight, py ^ Z, is satisfied. 

The other condition concerns the separation of the singularities 7^ fy for i 7^ j. The four singu¬ 
larities are governed by a number of parameters: the discriminant X> 


(44) 


:= (zi -hZ2Z3)(Z2-hZiZ3)(Z3-HZiZ2)(l-HZiZ2Z3), 


the additional auxiliary variables 


(45) 


T := 1 -h vlv 2 — 


' + ^2)^3 


V := 


4 A 


(1-b 371)2(1-M;2)2(l-h 1^3)2 


:= 


( 1 - 


— 4 l 32 r 2 




.2'i2 


SO that 

(46) 

A^ := (1 131132 - V1V2, - 372133) (1 - 131132 - 131373 -h 372133) (1 - 371132 -h Z 3 i 373 - 372373 ) (1 -h V1V2 + V^V^, -h I32Z33), 


and 

(47) 

A 2 := (1 - V1V2 -h 131373 -h 132133) (1 -h 331132 -h 33il33 - 332333) (1 -h 37i332 - 33i373 -h 372373) (1 - 33 i 372 - 37i333 - 333333). 


The singularities have the explicit forms 


^ T-hA ^ T-A ^ T-hA ^ T-A 

(4^) 3-1 TTI' ^2 - 2 z,iz; 2 ( 1 -333)2' ^4 - 


2331372(1 - 333)2 


2331372(1 -h 333)2 


2331332(1 -h 333)2 


Furthermore we note the relations between the singularities 


(49) 


^1^2^3^4 = 1, C3=zICi=C2^, ^4=zi?2 = ^r'. 


The sub-resultants of the numerator and denominator polynomials in are proportional to 
and 33i332333r. The Curie point 1I28I1 . the three equivalent variations of the Neel point and 
the degenerate case of T = 0 account for the vanishing of the first of these sub-resultants, whereby 
a single cancellation of a common factor in the numerator and denominator occurs and reduces the 
number of singularities by two. Examining the sub-resultants of the numerator polynomial and its de¬ 
rivative we find 333332(1 — 333 ) 2 a 2 and V1V2P — 333)^. Proceeding in the same way with the denominator 
polynomial and its derivative we find v\V2l^ and 333332. These last two cases include the Disorder point 
| [3d) l and its three equivalent versions, as well as the case of T = 00. Thus in contrast to the generic 
case the three critical points and the boundary points T = 0 , T = 00, correspond to situations where 
the singularities coalesce in a pair-wise marmer in the following way: 








Curie or Neel Point 

(50) = 1> 1^2! = K3| = 1> 1^41= zi 

23 

or 

(51) |^2| = 4>l^ll = K4|=l>K3l=zi 

23 

T = 0 

(52) = ^3 ^ ^2 = ^4 

Disorder Point 

( 53 ) Ki| = K2| = ^> 1 >K 3 | = K 4 | = |z 3 | 

T = 00 

(54) ft = ^3 ^ ^2 = U 


The reason why there are only three independent variables is that the integration contour defining 
w„ from the weight can be contracted or dilated from |^| = 1 to one of the nearby singularity mod¬ 
uli (keeping the fixed singularities at ^ = 0, oo unchanged) which has the effect of normalising the 
remaining three by this particular one. □ 

Our first result of the identification of a Gamier system with the spin-spin correlations is a simple 
recurrence relation for the Toeplitz matrix elements. 


Corollary 2 ([^). The Fourier coefficients or Toeplitz matrix elements Wn satisfy the fourth order, linear 
homogeneous difference equation in the index n 


Wn -2 


(55) {n - 3 )vlvl{l - vl)^w„_3 - 2V1V2T V3+ {n-2){l + vl) 

(« — 1) + 4:V2vj + ^3 “ 2(n — 1) + V2 — 6v\v 2 + v\v2 + 1^3 


-|-(m - 1) (1+ 4 v^V 2 + -h v^) 


Wn-1 


— 2 v-iT’ 2 r 


V3 + n(l + V3) ZVn + {n + - V^yiVn+i = 0 . 


In addition for n = — 1,3 the order drops to third order. One could take any contiguous set of four elements 
including wq as the initial values and iterate in either direction. 

Proof. This follows from an adaptation of Eq. (4.36) of m to the case where £4 7^ 0 


(56) 


E (“1)^ [(P “ 0^3-/ - "23-;] ^p-i = 0/ P e Z, 

;=-i 


and the spectral data (eg = 1) 

,(l + p§) {{vl + vD vl-vlvl- 1 ) 

ViV 2 il-vl)^ 


(57) 

£1 = £3 

= 2- 

(58) 

£2 

("t 

(59) 

£4 

= 1, 

(60) 

mo 

= 0, 

(61) 

nil = ^3 


(62) 

m2 

- 8 ’ 


{v\ + 41 ^ 2 "! + "2) "3 “ 2 {v\v\ + V2v\ — 6 v\v2 + v\ + v\ + v\v2 “H 4 v\v\ + 1 


"p2(l -"3) 


2t2 


"3 + " 3 -"l" 2 ~l) 

z;i 1^2(1 - 

"3(1+ "3) 

(1-1^2)2 • 
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Note that wq cannot be determined by this recurrence by any set of zv,, either n < 0 or n > 0 and is an 
arbitrary normalisation of the weight. □ 


However the theory of discrete Gamier systems developed in m furnishes a direct recurrence 
relation system for the Toeplitz determinants themselves, and our second and more significant result 
is such a system for our application. Prior to stating these we need define appropriate co-ordinates by 
choosing one of fhe singularities, say ^4 without any loss of generality 


(63) 


54 0 «K 4 )' 


Sn — f i — 1/ 2, 3 . 


The variables &n{z) and coh arise in the parameterisation of fhe derivatives of the bi-orthogonal poly¬ 
nomials and associated functions with respect to z, and is termed the spectral structure of the isomon- 
odromic system. Further details can be found in j^, l^u] . 

As a preliminary sfep we need to make the following definitions in order to render the results in 
the simplest possible form. 


Definition 1. Let us define four auxiliary variables 


(64) 

7 ^l 

(65) 

7^2 

( 66 ) 
7^3 

(67) 

'JZ4. 


(F-A) 1 

2 z;iU2(1 -Fu 3 ) 2 ^” 

_ + pi + 

2 z;iU2(1 

_ + pi + 

2viV2{1 - 

(F-A) 1 

2 z;iU2(1 


( 21^1 172 ( 1 +^’3)^) 

( r+A y 2 

\ 2 viV 2{1 + V 3 )^ J 

( r + A 2 

V 2171172(1 


(2171172(1-173)2) 


f F-A y 3 / F-A y 

^ 2171172(1 + 1 ^ 3)^7 ^”"^125^11^2(1 + 1^3)^7 

( r+A y 3 +( r+A y 

V 2 Z 7 ll 72 (l +173)27 V 2 l 7 ll 72 (l +173)^7 

( r+A y 3 +( r+A y 

V2z7ii72(l-153)^7 V2 i5i152(1-173)27 

/ F-A y 3 / F-A y 

V 2 z 7 i 172(1 -153)^7 V 2 i 5 i 152 ( 1 -153)^7 


Furthermore we define another set of four auxiliary variables 


(68) 5 i := ( 2 Fi 73 + A(l + z 7 i)) 


(1-153)^ 
.(1 + 153)^ 


(F + A)-(F-A)/; 


+ (2Fi73-A(l + i72)) 




(69) 52:= (2Fi73 + A(l+z7i)) 


(1-153)^ 
(1 + 153)^ 


(F + A) 2 -(F-A) 2/2 


(2FZ73 - A(1 + 172)) ^ ^^ 2^1 _ |i-p|^(F - A)2/; 
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( 70 ) Ss := (2ri;3 + A(1 + 




(r(3i;i - 2 z;3 + 3) + A(1 + v^)^) (F + A) 
- (r(3z;i + 2 i;3 + 3) - A(1 - v^)^) (F - A)/^ 

+ (2Fz;3 - A(1 + vD) (F(3z;i + 2^3 + 3) + A(1 - v^)^) (F + A)fl 




(F(3i;i - 21^3 + 3) - A(1 + 1^3)2) (F - A)/3 


(71) 54 :— ^2Fi?3 + A (1 + ^3)^ (f(3i;3 + 2z;3 + 3 ) — A (1 — 173)^^ (F + A) 

- (F( 3 i^i - 2 z ;3 + 3 ) + A (1 + vs)^) (F - A)^f; 

+ (iTv^ - A (1 + z; 2 )) (f( 3 i ;2 - 21^3 + 3 ) - A (1 + v^)^) (F + A)V, 


(1 +^’3)'^ 


(f(3i;2 + 2z;3 + 3) + A(1 - v^)^) (F - A) V,? 


Then one has the following three-variable generalisation of the additive member of the Sakai 
classification m otherwise known as the "fifth discrete Painleve equation" . Subsequently we will 
give explicit relationships between the co-ordinates introduced above and the Toeplitz determinants. 

Corollary 3 (|j^). Assume that the singularities pairwise distinct i ^ j, i.e. that A 0 , A 7^ 0 

and that 81,62 ^ 0 . The set of variables {fl,/„> fn>Sn>Sn>8n}n=o satisfy the system of coupled, first order 
non-linear difference equations in n > 0 . The first set of three constitute the first members of the Lax pair 


(72) i! , ?L. /n/nVi = 


[7^l - n] 


TZ'i — n -\ —2~2 


2P3A 


.^ (A^-4ptP2P3(l + P3) ) _ A 
Z;^i;^(l + I^ 3 )' I r + A 


(l + ^^ 3 )^' 


(73) 


(1 ~ ^3) (r -|- A) 2 r2 ^ 
(1 + z;3 )2(F- 


[R-i — n 
[ 7 ^ 2 -n 


72-4 — n -\ —j—j 


2mA 


^1^2(1 -^3) 


(A^-S 4 z;f 322^3(1-i;3)^) 
F + A 


-A 


72-2 — n- j 


2V3A 


(A -4i;ii;2i;3(l + ^3) ) , a 


viv'iii + vs)^ r - A 


[724 - n] 


724 — n —j—j 


2P3A 


,^ {A^ + Avp^V3{l-V3)^) 
vivi{l-V3)^\ r + A 


2^ ^ ' 
-A 


(74) 


r + A ^3 ^3 

F-A 


[723 - n] 


m+i = 


723 - n - -j 


2P3A 


vivi{i-v 3 r 


(A^ + 4 p^u^P 3 ^( 1 -i; 3 )^) ^ 


[724 - n] 


724 — n -\ —j—j 


2P3A 


^ (A +4pt 132^3(1- 133 ) ) 
vivi{l-V3)'^\ r + A 


The remaining set of coupled, first order non-linear difference equations constitute the three members of the 
second Lax pair and are given by n>l 


(75) gl+gl-i- 


2viV2{'l - V3)^ 


2 n-l + { 2 n- 3) 


(1 + 5^3)^ 


n +1 


52 


2z;ii; 2(1 - 5i 2z;ii72(l - ^2 


= 0 , 
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( 76 ) 


, , 81 ^ 3 ( 1 + ^i) 


^ 1 + 4v^V2 + vfvj — 2(vf + 'v\ 


6v\v\ + v\v\ + v\v^v^ + iv\ + 4vfv2 + 

^ 

_ n + 1 _^ n S3 ^ 

+ 323 ) 2(1 - 1^3)4 5 i (1 + 133)252 


(77) gl+gl-1- 


2viV2{1 - 133)2 


2 n - 1 + ( 2 m - 3 ) 


(l + ^’s)^ 

M + 1 


2 z 3 i 132(1 - 13^)2 5 i 

These recurrences are subject to the initial values for j = 1 , 2,3 


- 2 mi 3 iZ 32(1 - 323)^^ = 0. 

<62 


(78) = 


1 _ (1 - ^ 3 ? 


(1 + ^3)^ 


(r-A) ^ (r-A) (r - A - 


vjv^(T+vfj^^ 2 + 4 i,++i+j,g )4 


(1-333)' 


1 if.i - 




(r-A) 


DiD 2 ( 1 -I! 3 ) 


jW—2 “I” 


(r-A) 

ivlvKl-v^)^^ 


r-A-ltgjilr 


3X3_i 


lv?T 

VlV2{l-vl)^ 


Wq + Wi 


(79) /o 


(i-^3)2(r + A) 

(1 + Z33)2(r-A) 

(r+A) (r+A) 

Ul 112(1 + 333)^ ^“2 + 4 i;Ji;^(l + i; 3)4 

X -p- 

(r-A) (r-A) 

uii; 2 ( 1 - 333)^®“2 + 4 uJ»^( 1 -z; 3 )'‘ 


r + A 


4(1-333+333) r 
(1-333)^ 


lV_i 


r-A- 


4 ( 1 -P 3 +d§) 

(1+333)^ 


If-l 


z,iV2(l-vi)2^0 + ^1 

v,v2Ti-vI)2^‘’0 + “"1 


(80) fo^ 


and 

(81) 

(82) 

(83) 


(r + A) [131332(1-133) 


(r+A) (r+A) 

-2^11-2 + ^; 5;;3 - 


-fl 32 (l-i 33)4 


(r + A - + -■», 


(r-A) 


(r-A) 

Ul 1)2(1-333)' 


rro_2 + 


(r-A) 

4 d(p|( 1 -D 3)4 


r_A_ 4 (i_^^r 

(1+333)^ 




1 ^ o Par _ ^ 

° 131172(1-133)2 Ifo' 

2 ^ 133(1 + 13§) I73r aiQ (1 + 173 + I3^)r wi ail W2 

° (1 — 133)2 13iZ 32(1 — Z32)2 a;_p I3iZ32(l — 13^)2 ZMq '^-1 

3 ^ 2 ^ - ^ - 2 (^~^3 + ^ 3 )r 

K3-1 M30 131132 ( 1 - 133 )^' 


Proof. We will first permute and re-label the singularities in the following way ^i = ^4, ^2 = ^3/ 
^3 = ?!/ ?4 = ^2 so that we can rescale the contour radius in order to normalise with respect to ^4. 
From an adaptation of the first result in Proposition 4.3 of [ gn] the first Lax pair can be written for 
; = 1 , 2,3 as 


j j [^jgl + ^jgl + ^]gl + ^j{'^ + mo)-ne4 

?/■ {gl, + m) + [gl - m) + {gl + Mil) + ?^ - Me 4 

[?4 {gl + m3) + {gl - m2) + {gl + mi) + - Me4] ' 

Substituting the evaluations of the singular points and employing the definitions |[6^-l[6^ we 
simplify the resulting expressions and arrive at 
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Proceeding we now define [T = {a,b ,...>,jgT f i„ as the n-th elementary sym¬ 
metric function in the variables the Vandermode determinant A[T = {a,b,...}] = 

nr<se{fl,i7,...}(?s “ ^r)- Adapting now the second result in Proposition 4.3 of | j^ the second mem¬ 
ber of the Lax pair is given for j = 1 , 2,3 by 


(85) gL+g’n-l + (-l)^(w - l)e 4 -/ - 

.(,y+i( ,1 , S4-y[123]A[123] - 54-y[234]A[234]/i + 54-;[134]A[134]/„^ - 54-y[124]A[124]/3 

A[123]-A[234]/1 + A[134]/2-A[124]/3 


(- 1)^4 


S3_,[123]A[123] - S3_/[234]A[234]/1 + S3_/[134]A[134]/2 - S3_/[124]A[124]/3 


= 0 . 


S3 [ 123 ] A[ 123 ] - S3[234]A[234]/i + S 3 [ 134 ]A[ 134]/2 - S 3 [ 124 ]A[ 124]/3 
Using the definitions llMll -([7i|l and the evaluations of the other data we deduce 

The first set of initial values, ([78)- ll8ot , are ratios of the spectral coefficient (see Eqs. (2.6) and (2.10) 
of 133 ) 

(86) 0„(O = 1 

Yn-\-l\^) ^n+1 ^n+1 ^n+1 

appropriately specialised using the data |[37ll- (l62ll and evaluated at ^ = fy for j = 1 ,..., 4 . The second 
set, l|8i]l-l|83, can be derived from a third spectral polynomial U{Q using the generic relation 

(87) no(0 + v(0 = ^ (1 + [2V{0 - K^uiO 

and the fact that = [^^](no + which follows from the parameterisation of n„. The polynomial 
U has the explicit form (which is an adaptation of Eqs. (4.37)-(4.39) of m for £4 0) with coefficients 

(88) uq = le^wi -h nt^WQ, 1/3 = niQZVo, 

(89) M; = (-l)^"'m3_;w;o+ 2 ^ (-1)^“''[(/- r)e3_r - m3_,] W;.,, 1 = 1,2, 

r=-l 

in terms of a contiguous set of initial moments, essentially the same set that defines the initial values 
of the recurrence in Corollary Thus U encodes these initial values in an alternative way. □ 


Under the conditions applying to Corollary [3] the primary variables can be recovered from the 
ifLgit} variables using the formulae given in the following corollary and is a direct consequence of 
the theory in |9|, |4T] . 

Corollary 4 (||^, [^). Assume the conditions of Corollary^ The Toeplitz determinants < [31) 1 and l [^ , abbre¬ 
viated as {o'o,oO'm,n)^ —■ 1n/ computed using the recurrence relation equation 

^n+lhi—l 


(90) 




= 1 - r„r„, n>l, 


with the initial conditions Jq = 1/ li = ^o- Iw turn the pair of reflection coefficients r„,r„ are computed firstly 
using the recurrence relation for r„ 

r„+i _ 2m;iz;2(l - Si 


rn n + 1 S 2 ' 

The coefficient fn is computed via another coefficient A„ which satisfies 


n > 1. 


(92) 


(n -|- l)A„+i — nAn — ~ 2 w _ _ - j— g„-\- ■ 


n + 1 


S3 


viV 2 {l — v^fi 2 z ; ii ; 2(1 — 1 ^ 3 )^‘^ 1 ' 

and the generic relation — A„ = Alternatively one can use 


(93) {n + l)A„+i-nA„ = gl- 


2viV 2{1 - vsfi 


2n 




jn + l) S 2 
(1-h 1^3)2J ' 2 viV2{1 - V3fi Si' 


-H 


along with A„ — A„_i = r„r„_i to determine r„. The initial conditions are ro = ro = 1 and Aq = Aq = 0 . 
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Proof. From the general theory the spectral coefficients can be parameterised, in the case of four finite 
singularities, in the following forms (see Eqs. ( 2 . 6 ) and ( 2 . 10 ) of m) 

fn{0) a ^ ^ + 1 + 


(94) 


&n = -ne^- 


^n+1 


^n+1 


^n+1 


<^n+l( 0 ) 

and (see Eqs. ( 2 . 7 ) and ( 2 . 11 ) of ED) 

(95) = -ne4Z“^ + + {col “ + {eol + + (1 + \mQ)z^. 

The parameters introduced above can, in turn, be related to the leading or trailing coefficients (with 
respect to z) of the bi-orthogonal polynomials via the formulae 


(96) 

= ne^-t 


rn 

(97) 

K = -{n + l)ei - 

(98) 

col = ne3 - 

(99) 

col = “^1 -fni + 


{n + — (n — 1) ( A„_i + 


^n—l 


rn+2 

^n+1 


-A. 


n+2 


Tn 

+ (n + mo)A„, 


(n + l)A„+i — n I A„ + 


hi+i J \ 

^n+2 


rn+2 

rn+l\ 


However these parameters are also given by the dynamical co-ordinates by the formulae 

rn _ M + 1 + mo S 3 [123] A [123] - S 3 [234] A[234]+ S3[134] A[134]/2 - S3[124]A[124]/3 


( 100 ) 

( 101 ) 
and 

( 102 ) 


^17 + 1 


^1 = {n + l + mo) 


K = -{ri -t -1 -I- mo) 


A[123] - A[234]/i + A[134]/2 - A[124]/3 

S 2 [123] A [123] - S2[234]A[234]/1 + S2[134] A[134]/2 - S2[124]A[124]/3 
A[123] - A[234]/i + A[134]/2 - A[124]/3 

Si [123] A [123] - Si[234]A[234]/i + si[134]A[134]/2 - si [124]A[124]/3 


A[123] - A[234]/i + A[134]/2 - A[124]/3 

Combining and ( [ 9 ^ we note 

(103) col-dl = nei + {n + 1 + mo)K+i - {n + mo)A„. 

The first of these, dioob . gives the formula ( [gi) . From the preceding relation i |i 03 | i and (1 102 b we deduce 
the formula Ijg^. There are also other ways to determine fn, such as using ( [ 98 ) with ilioob , and this 
yields ([ 93 ^. In addition to the forgoing specific relations we have the general identities 


(104) 


^n+1 — ^n+l^n/ '^n+1 


which allow us to close the system. Finally the Toeplitz determinants satisfy the second order differ¬ 
ence equation 


( 105 ) 

which is also a general identity. 


Ifi-l-lln—1 

J2 


= 1 - r„r„. 


□ 


In the remaining part of our study we discuss two important special cases that arise from the 
general formulae above by taking appropriate limits as given in ([ 37 I 1 and 
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2. Correlations along a Row/Column of the Rectangular Ising Model. The row and column correla¬ 
tions of the rectangular Ising model also have a Toeplitz determinant form (see ¥111.2.28,29,30 of lISTl '). 
Because the row and column correlations are given by each other under the exchange Zi zi we will 
treat only the column correlations without any loss of generality. These latter correlations are given by 

(106) (c^0,0C^0,n)'^ = det(n;y_n.)o<;,J:<N-l/ 


where the Fourier coefficients Wn 

d6 (1 -F ixia2 — («! -F ^2) cos(F)) cos(n0) — (ai — X2) sm(0) sm(n0) 

l—n 2/1 


(107) 


= £ 


are defined by the weight 

(108) w { C ) = 
Two new co-ordinates are defined by 

(109) 


— 2 oci cos( 0 ) -F af) (1 — 2 a.2 cos{6) + a|) 


( 1-^20 J 


1 — Zi 1 1 — Zi 

= Z2:r--/ Ol-l = 


■ 1 -F Zi Z2 1 -F Zi 

The parameters aj, oc2 are also related to the set S, C and S, C by 


(no) 


«1 = %^(C-S), a2 = ^(C-S). 


which is entirely analogous to the definitions of McCoy and Wu but differs because of conventions for 
the lattice co-ordinate system. The modulus is related to these parameters by 


(ill) 


1 — OL1OL2 

a.2 — 


For 0 < k < CO the parameters satisfy the inequalities 

(112) 0 < ai < ^2 < 1 < 0-2^ < A: > 1 , T < Tq, 

(113) 0 < «! < ^2 ^ < 1 < ^2 < k <\, T > Tq. 

Lemma 1. The Fourier coejficients satisfy the symmetry 

(114) W„{Ul,U2) = 'W-„{U2,0il) 


As indicated by l|37f the result (I106II and {loj) can be found from the diagonal correlations on the 
triangular lattice by setting K2 0 followed by the relabelling ^ K2. This is still within the M = 4 
semi-classical class of the general case even though the number of independent variables has been 
reduced by one. The singular points are now 

(115) (^i/^ 2 ,^ 3 /C 4 ) = (ai"£a2,a2£ai). 

Even though this is merely a specialisation within the original system we will record the full details of 
the final recurrence relations for the convenience of the reader since they are new and of separate and 
intrinsic interest. 


Corollary 5. The Toeplitz matrix elements w„for the column correlations of the anisotropic square lattice Ising 
model satisfy the fourth order, linear homogeneous difference equation in the index n e Z 


(116) 2 aiiX2{n — 3)a;„_3 — (1 -F OC1U2) [{ 2 n — 5 )ai -F ( 2 n — 3)^2] zVn-2 

(n — 2 )a.l + na.2 -F (n — 1)(1 -F aiA;2)^ 


W„-i 


— (1 -F 0.1X2) [{ 2 n — \)xi -F ( 2 n -F 1)^2] + 2 xiX2{n -F l)zv„^i = 0 . 

One cotdd take any contiguous set of four elements including ivq as the initial values and iterate in either 
direction. 
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And the corresponding specialisation of the coupled non-linear recurrences for the Toeplitz deter¬ 
minants is given by the following result. 

Corollary 6. Assume that (1 — a^) [l — cc^fn] — (1 — [fl — oi^fn] 0 and ai(l — a^) [l — oc^fn] ~ 
0:2(1 — 0.2) [f^ — ^ 0 . The set of variables {fl, fn> fn> sh' si' sV\n=Q column correlations of the 

anisotropic square lattice Ising model satisfy a system of coupled, first order non-linear difference equations in 
n > 0 . The first set of three constitute the first members of the Lax pair 


( 117 ) 




~ + “iS'n ” 2 + 2^^! + 2a\ 


[ml + “kn + - n] 


ml + ‘’^Isn + — W + 2 + 2«2 “ + “l)(l “ “ 2 ) 


( 118 ) 4fnfLl = 


[4sl + « 2 gn + ml - ^4 +1] 

[ml + + 4gl + ^^2 - "] 


4.gl + ^Igl + ml - (« + \)4 + 2 - ^(1 + “ 1 )(1 - “ 2 ) 

ml + 4 gl + -n+l + ^xj- 1^(1 -h af)(l - 4 ) 


(119) : 4 fnfn+l = 


X2 


[4gl ‘^ign ~b‘^Ign ~ ^4 ‘^igl ‘^igl 2 )^ 1 + 2 + 


di 




[ml + 4gl +‘^igl +‘^2 ”] x2gi-\-4gl + 4gl~^ + i + i‘^ 2 ~^i^ + 4 )(^~‘^ 2 ) 


The remaining set of coupled, first order non-linear difference equations constitute the three members of the 
second Lax pair and are given by 


( 120 ) gl-Lgl_i 


il±^[(2n-3K + (2n-l)a2] 

1 / [l-«:2/n] -“2(l-a^) [/n-at/«]] 

^ [( 1 -« 2 )[i^ 2 ^ 2 ](i^ 2 )[^ 1 ^ 2 ^ 3 ]] 

(1 -a^) [(1 -hai(ai -^0:2)) - (ai -h 0:2(1 + af))a^/^] 

- (1 -“ 2 ) [(l + «2(ai +a2))/n - (a2 + ai(l + a2))'^i/«] 

[«l(l - 4) [1 - ^t/n] - ^ 2(1 - 4) [fl - 4fn]] 


( 121 ) 


?l + gl-1 


fto — OC^ 


-\-{n-l) - o:f-h 0:2-h ( 1 -I-0:10:2)' 


OC1CC2 OC1OC2 L 

(1 - x\) [(0:1 -h 0:2(1 -h 0:^)) - (1 + o:i(o:i -h o: 2 ))o: 3 / 2 ] 

- [(“2 + o:i(l -Fo:|))/l - (1 -ho:2(o:i -Fo; 2 ))o: 3 / 3 ] 


[o;i(l - a.\) [1 - o:f/ 2 ] - 0:2(1 - 0:^) [/l - xjf?,]] 

(1 - x\) [(1 -h o:i(o:i -h 0:2)) - (0:1 -h 0:2(1 -h x\))xlfl] 


n -1- 1 


- (1-0:2) |‘(i .^^2(0:1 -F 0:2))/! - (0:2 0:1(1 -ho: 2 ))o: 3 / 3 ] 


0:10:2 


[(1 - 0:2) [1 _ ^2y2] - (1 - 0.1) [/I - 0:2/3] j 


= 0 , 
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(122) gl + gl_i 


(1 + a:i«2) 

2DLiCi2 


Kl0i2n 


[{2n — 3 )ai + {2n — l)(x.2] 

[(1 - al) [1 - alf^] - (1 - Oil) [fl - ] 

[ai(l - ot\) [1 - oc\fl] - a2(l - 4) [fn - 4fn]] 

(1 - 4) [(ai + 0i2{l + 4 )) - (1 + + ‘>^ 2 )) 4 fn] 

- (1 - 4 ) [(«2 + ai(l + 4 ))fn - (1 + ‘^2(^1 + 0 i 2 )) 4 fn] 


KlOi 2 [(1 - 4) [1 - 4fi] - (1 - 4) [fh - 4fn] ] 

These recurrences are subject to the initial values for j = 1 , 2,3 
(123) 

j Oil 44 ^ 2^-2 — *i('^ 2 (l + 4) + <^1(3 + 4))^-i ~ (“2(1 ~ 4) ~ 4))'^0 + 2 aia 2 ifi 


= 0 . 


/o = 


0C2 4 :ai 4^-2 ~ '^2(®2(1 + 3 aj) + ai (3 — i4))'w_i + (ai(l — <4) — ^2(1 ~ ®i))wo + 2aia2ifi' 

2, 


(124) /o = 

(125) 4 

and 


,2 1 4aia2if_2 — («! + ^2 + aia2(3ai + a2))w_i + a2(^i “ '^2)(1 + aia2)roo + 2aia2it'i 


(126) 


(127) 

^0 = 

(128) 

^0 = 


4 4aia2W_2 — a2(3ai + 0 L 2 + aia2(3ai — ol2))w_i + (ai — 0C2){1 + ciia.2)wQ + 2aia:2Wi' 

1 4a2a:iiV-2 — {3oi\ — 1 X 2 + aia2(3ai + a2))ro_i + ai(ai — 0 : 2 ) (1 + aioC2)iVo + 2a^a2Zi7i 
4oi2 4ai4^-2 — a2(3ai + ^2 + aia2(3ai — a2))zi^_i + (ai — < 12 ) (1 + aia2)wo + 2oiiOi2Wi' 

(ai(l — (4) — ^ 2(1 “ + 2aia2roi 


^0 = 


'2.0C\DL2'Wf) 


2{4 ~ ^2)W-iWo + (ai + 30:2) (1 + aia 2 )w_iwi 

— 4 aia 2 if iif 2 + ('^i ~ ®2)(1 + oiiOL2)v4 + 2aia2ii’oifi 


2aia2if_iZ(;o 

4aia2it'_2Zt’o “ 2oli(X2w4 — (ai(3 + a^) + a2(l + 3aj))zi;o“’-i 
2aia2if_iZ(;o 

The primary variables can be recovered from the the {fli, 4 n} variables using the formulae given in 
the following corollary and is a direct consequence of the theory in 0. Ell- 

Corollary 7. The Toeplitz determinants l lio6b and {107) , abbreviated as (uo^o^o,n)'^ =• ^N/ computed using 
the generic recurrence relations of Cor. g] with the initial conditions Iq = 1 , Ii = wq. In turn the pair of 
reflection coefficients r„,r„ are computed firstly using the recurrence relation 

r„+l _ 0ii0C2n (1 - 4 ) [1 - > 44 ] - (1 - 4 ) [fn - 4 fn] 


In n + 1 ai{l - 4 ) - 44 ] - <^24 - 4 ) [fn - 4 fn] 

Secondly ?„ is computed through another coefficient An ivhich satisfies 


n > 1. 


(130) (n + l)A„+i-nA„ - 


3 (ai+a2)(l + aia2) 


0.10-2 


1 


(1 - 4 ) [(“1 +0.2(1+ 4 )) - (1 + +'^2))a2/^] 

- (1 - ol) [{02 + ai(l + ol))ffi - (1 + a2(ai + a2))ai/n] 


“i“2 (1 - ol) [1 - c4ffi] - (1 - ol) [ffi - olffi] 

and the generic relation A„+i — An = in+il’n- Alternatively one can use 

(131) (n + l)A„+i - nXn = gn - [{ 2 n - l)ai + ( 2 n + 1)^2] 


20i02 


+ 


n + 1041- ffif) [1 - ojffi] -02(1- ol) [fn - 4fn] 

0102 (1 - ol) [1 - olffi] - (1 - ol) [/I - olffi] 


along with the previous generic relations. 


17 



























3- Diagonal Correlations of the Square Lattice Model. The diagonal correlations of the rectangular 
Ising model also have a Toeplitz determinant form (see Eqs. (3.29,30,31) of t2i1 ~) given by 


(132) 

where the Fourier coefficients iVn 

(133) IVn = 
are defined by the weight 

(134) ^{0 = 


(c^o,oc^n,n) — det(z(;y_j-)o<;,/c<N-i/ 

d 9 kcos{n 9 ) — cos{{n + 1 ) 0 ) 
-Tz 2 n — lkcos{ 9 ) + k^) 


L 


(i-fc-'O J ' 

and where k = SS, S = sinh 2 fCi,S = sinh2fC2. This arises from ([3^ as Z3 —?■ 0 or ^3 —> 1 and thus 
a = 4ziZ2, b = 0 ,c = (1 — Zj) (1 — Zj), so that k = a/c. We will also use 


(135) 


^ ^ 1 ^ (l-zf)(l-z^) 
k 4 ziZ 2 


Under this limiting process the singularities behave like 


(136) 

(137) 

(138) 

(139) 


Cl = -{V3-1) +0((U3-1) ), 

^2 = a + 0((U3-l)l), 

^3 = l+0((i;3-l)'), 

U = ^(^3-l)^ + 0((i;3-l)3), 


and so two of these merge with those at 0 ,00 leaving two at finite, non-zero locations. Therefore the 
three variable Gamier system reduces to the single variable Painleve VI system. We first consider the 
limit of Corollary |2l 

Corollary 8 ([H). Under the limiting process Z3 —>■ 0 or ^3 —>■ 1 Eq. 155)1 /or the Toeplitz matrix elements Wn 
reduces to the second order linear difference equation 

(140) a( 2 n — 3)io„_2 — 2 (a^-|-l)n — 1 io„_i-|-a( 2 n-|-l)zo„ = 0 . 

This is identical to Eq. (2.18) of tii. 

Next we check the limit of the nonlinear system in Corollary |3] 

Proposition 1. Through the limiting process Z3 —?■ 0 or 03 ^ 1 the dynamical variables 

(141) fl = {v3-lM, fl = Tl /3 = (u3-1)-"F„/ 

(142) gl = {V3 - gl = {V3 - gl = {V3 - 1 )-^gI, 

where F^, Gn = 0 ( 1 ) in this limiting process. Let us define the change of variables 

n2 _ /„~ 4 r /^2 


behave as 


(143) 


= a f„, G„ = -gn. 


Under the conditions gn + na ^ 0 and g„ + {n+ j)oc ^ — jOc 0 the dynamical equations reduce 

under the limit to the pair of coupled first order difference equations 

[gn +{n + l)a.-a~'^] g„ + {n + - ^a-'^ 


(144) 


® fnfn+l 


(145) gn+gn-i+2noi ^-2 [ol + ol '^) + {n + ^) 
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[gn+noc-f gn + {n + - jOi 

1 (a2-i) 1 


1 -/n 


+ (” + 2) 


. a [of — 1) 


cU- fn 


= 0 . 

















This is precisely the discrete Painleve system in the Sakai scheme otherwise known as "discrete Painleve 
V". These are subject to the initial conditions 


( 146 ) 


/o = a 


W-i + OCWq 
aw -I + Wo' 


1 / w_ 


So - 2 


Wo 


Wo W-1 


Proof. We note some preliminary expansions of a number of auxiliary quantities appearing in the Lax 
pairs ([7^-([74]l and 


(147) 

( 148 ) 

(149) 

( 150 ) 


2Au3 

2{A'^ - 4 vlvlv 3 {l + Vs)^) 

v\v\{l + V3)<= 

T + A 

2Au3 

2 (A^ - 4 i; 2 u 2 i; 3 (l + 1 ^ 3 ) 2 ) 

v\v\{l + V3)^ 

T-A 

2Av3 

2 (a2 + 4 i;^z;^iJ 3 (l - 1^3)2) 

vlv\{l - 03)6 

T-A 

2^V3 

2 {A^ + 4 v\vlv 3 {l - 03)^) 


v\v\{l - 
along with 

( 151 ) 

(152) 

(153) 

(154) 

and 


r + A 


-A 


+ A 


■A 


-A 




2(1 -a2) 


128 


(U3-I) 2 + 0((i;3-1) ^), 


= --^(u3-l)-« + O((u3-l)-0, 


=-2(1-a^)(u3 - 1) ^+0((U3-1) 


= 4 ^n (^3 - 1 )^ + 0 ((U 3 - 1 ) 3 ), 


7^2 = 


« ^sl + ^ ^8l + <^ ^ 


■ 0 ((U 3 - 1 )), 


7^3 = ^(^3-l)-' + 0((z;3-l)-a 


7?.4 = 


+ ^^Sl + 


+ 0 ((U 3 - 1 )), 


(155) 

( 156 ) 

(157) 

(158) 


Si= 2 {l-v\f{l-vlf (z;3 -1)2 + 0 ((u3 -1)3), 

^2 = -4(1 - 1;2)3(1 _ ^2)3(4 _ _ 4)2 ^ _ 4)3)^ 

^3 = 16(1 -u?) 3 (l-z ;^)3 I - of)fl- offl (i, 3 _ 1)2 + 0 ((U 3 - 1 )'‘), 

54 = 8(1 - [1 - (1 - - aVnj (^3 “ 1)"^ + 0((z;3 - 1)3). 


Using the scalings i |i42| i in ([7^, and assuming + (xG^ + oc^G^ ^ 0 and aG^ + ofG^ + a^G^ — 2(1 — 
a^) 7^ 0 we deduce from the leading order in this expansion that 

(159) (aG^i — 4 n)(aGl — 4 n — 2 ) = 0 . 

Treating in a similar manner we conclude, again under the above conditions, that 

(160) (iaG 3 + l)(iaG 3 + l) =0. 

Which of these two sets of solutions are relevant can be settled by the initial conditions, which we 
undertake subsequently. The other member of the first set of the Lax pair, 1 (7^ , now has at leading 
order the expression 


( 161 ) 


1 

a® 


6p2p2 _ 

^ ^n-^n+1 


x^Gl + ocGl + Gl] [oc^Gl + oc^Gl + (xGI - 2(1 - of)] \ 
[Gl + ocGl + a^G^] [ixGl + a^Gf, + a^G^ - 2(1 - off] ] ' 


19 



































Undertaking the expansion of the initial values given in |[78^- ll8ot and ll8ill -l(8^ we compute the leading 
terms as 


(162) 

(163) 

(164) 

(165) 

(166) 

(167) 
and thus 

(168) 

(169) 


/o = 


Wo 


4(aw;_i + Wo) 
w_i + azvo 
a^{aw-i + Wo) 

.3 16 w_i 


/o = 


/o =- 


a?{ciW-i + Wo) 


(V3-1)^ + 0((V3-1)^), 

+ 0(V3-1), 

(V3 - 1)~^ + 0((V3 - 1)~^), 


go = -(^3-1) ^ + 0((U3-1) 1), 


^0 = 


2aw-3Wo + Wn — 3 w-iW3 

_y_I 


(U3-I) ^ + 0((U3-l) ^), 


-I'l 


^aw_ 3 Wo 

-2 


^0 = --(^3-l) +0((U3-1) ), 


pi _ 
rn — 


Wo 


4 (ftw_i + Wq)' 


p^ — 
-un — 


W-l + aZVo 
a^(ciW-i + wo)' 


rj — _ 

if! — 


16 ro_ 


1 


a^{aw-i + Wo)' 


2 (10 


G^o = -, , 

a DC \ Wq W_3 


Wo 


Gg = —. 

" a 


From these initial values we see that the choice for the solutions to i|i39|i and lli6ob are 


(170) 


G,) = 


4 m + 2 


Gn = -: 


Now we turn our attention to the second set of the Lax pair and in particular expanding we find 
at the leading order that 


(171) 


1 1 2 ( 2 m- 1 ) 4 m 16 m 


G„ + G„_i - 


{l-a^)Fl 


a a 4(l-a2)Fi+a4F2-l‘ 


Equating this to zero and substituting the first of the solutions i |i7o| i for we compute F^ as 

/ , oi^F^ - 1 

“ 4 (a 2 -l)( 2 M + l)' 

Continuing we expand (j/Tjl to leading order and deduce the expression 


(173) 


4 m -2 


a^Fg — 1 


G„ + G„_i --+ 1 2n 2'ir3 i‘ 

a a a2F2 - - oc^)F,i - 1 


Utilising the second solution of i ji7o| i we solve for F^ as 

16 (a 2 F 2 - 1) 


(174) 


F^ = 


a2(a2 — 1)(2 m + 1)' 

Lastly we expand ( [76} and compute to leading order 

4 (l + fl: 2 )(M-l) 4 

a? a2 

4 (m + 1 ) 1 — a^F 2 4(1 —a 2 )M 


(175) G^ + G^_i + 


4F| - 1 




4(1 - a 2 )Fi + a 4 F 2 - 1 


= 0 . 


Finally employing the solutions ( |i7o| i and the relations i |i72| i and i liyqj i in the two remaining equalities, 
Ilr6ill and 1I175I1, we arrive at ^144} and 1I145I1 respectively. □ 
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Proposition 2. The diagonal correlations of the square lattice Ising model are governed by the generic relations 
of Corollary ^along with 


(176) 

and 


rn+l _ ^f„ 

rn 'i^-fn' 


1 — f 

(177) (n + 2 )A„+i - (n - i)A„ = -n(a + a"^) -gn + {n+ j) - 

OC OC tfi 

Proof Expanding the relation (j9^ as 1^3 ^ 1 and using the results of the previous proposition we find 
the leading order gives 


(178) 

Employing ||i72|| and 1I174I1 


^n+l 


— Ci- 


n + 1 
we deduce (176I1. 


ofPf, - 1 - - K^)Pl 

4(1 - a 2 )Fi + _ 1 


□ 


Remark 1. The Lax pair 1 I144I 1 and can also be verified directly from the theory in §4 of [4TI . From 
the weight 1 I134I 1 the configuration of singularities is given by 

0 a 00 

- 1/2 1/2 - 1/2 1/2 

from which we compute the spectral data as cq = = oi-+ = 1 arid mg = —l/2,mi = 

—a, m2 = — 1 / 2 . Employing the parameterisation of the spectral polynomials as 



(180) ^!l±i©„(^) = + 

Kfi 

( 181 ) n„{0 = l^^-{g„ + ^^)^-n-l 


and making the definition 
(182) 


a ^/n 


0n(a-i) 

©n(a) ' 


we can use Eq. (2.18) and Eq. (1.33) of m to deduce i jiqql i and i |i43| i respectively. This Lax pair differ 
from those of Prop. 4.1 of m only in the rescaling of the independent variable Alternative, but 
entirely equivalent, recurrence relations have been given in Corollary 5.3 of HI and Corollary 1 of [|4o|. 


Remark 2. The forgoing generic theory does not apply without modifications to the critical cases of the 
Curie, Neel and disordered points where A = 0 , A = 0 because of the undefined nature of the simul¬ 
taneous vanishing of numerators and denominators. There are two strategies one could adopt here: 
perform a local expansion about the critical points or adopt the more direct approach of evaluation of 
the primary variables from first principles. However we defer this task to a subsequent study because 
of the lengthy technical issues involved. 


Remark 3. We have not attempted to discuss the hypergeometric function evaluation of the Toeplitz 
matrix elements w„ in the general case of the various regimes: T < Tq, T^, or < T < Tq and T > 
Tq. Suffice it to say that one would expect the three-variable Lauricella D h5rpergeometric functions 
would appear in the diagonal correlations on the generic anisotropic triangular lattice whereas the 
third complete elliptic integral would describe the row /column correlations of the anisotropic square 
lattice, in addition to the first and second complete elliptic integrals. 
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